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Magnetic waves are a relevant component in the dynamics of the solar atmosphere. 
Their significance has increased because of their potential as a remote diagnostic 
tool and their presumed contribution to plasma heating processes. We discuss our 
current understanding on coronal heating by magnetic waves, based on recent ob¬ 
servational evidence and theoretical advances. The discussion starts with a selection 
of observational discoveries that have brought magnetic waves to the forefront of 
the coronal heating discussion. Then, our theoretical understanding on the nature 
and properties of the observed waves and the physical processes that have been 
proposed to explain observations are described. Particular attention is given to the 
sequence of processes that link observed wave characteristics with concealed energy 
transport, dissipation, and heat conversion. We conclude with a commentary on 
how the combination of theory and observations should help us understanding and 
quantifying magnetic wave heating of the solar atmosphere. 

Keywords: Sun, Magnetic Fields, Magnetohydrodynamics (MHD), Waves, 

Coronal Heating 


1. Introduction 


In spite of decades of advances in observations and theory, the solar coronal heating 
problem remains unsolved. A full understanding about the magnetic and plasma 
structuring of the corona, the nature of the plasma dynamics, and t he physica,! pro¬ 
cesses that create and maintain a hot atmosphere remain elusive (see Kunenjs, 


Withbroe fc Nove^ 19771: Kuperus et al. . 1981 : Narain &: UlmschneideiT 1996 : Aschwanden 


2005 ; Klimchud 2006t Hood , 2010t Parnell fc De Moortel , 2012 , for reviews ). 

The corona l he ating p r oblem originated more than seventy years ago, when 
Grotrian (193^ and lEdled ( 194,311 identified the presence of Fe IX and Ca XIV spec¬ 
tral lines in the light emitted by the corona, indicating the presence of fully i onised 
plasma at multi-mi llion degrees. I n a recent historical note, iPeter fc Dwivedil (120141) 
assign the merit to lAlfveii ( 194lh . who summarised six arguments to support that 
the corona is heated to such high temperatures. Since then, solar physics research 
has sought to identify the physical mechanism(s) that might b alance the thermal 
condu ction, radiation, and solar wind losses, first quantified bv lWithbroe fc Noves 
I I 977 I) . Researchers in the field agree upon the fact that the origin for such tem¬ 
peratures is of magnetic nature and that the energy source lies in the solar surface 
plasma motions. This available energy is then transported to the upper layers of the 
solar atmosphere where it gets dissipated. The exact physical processes and their 
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contribution to the heating o f the plasma remain largely unknown or unquantified 
( Parnell &: De Moortel . 2012 ). 

A number of mechanisms are believed to contribute to the heating: the di- 
rect dissipation of magnet ic energy by proces ses such as m agnetic reconnection 


([Stiirrock fc Uchidal. ll98lD . current cascades (lParkei| 1196,11). viscous t urbulence 

2004 ): the dis- 


van Ballegooiien . 19861) . or magnetic field braiding ( Peter et al. 


sipation of mag netic wave energy stored, t r ansported, and transferred to small 


scale s by waves ( Alfvenl . 1947 : lonson . 19781 : iHevvaerts fc Prie^ 198.ll : Goossenj . 


1991; Hood et ai 


19971) : or the mass flow cycle between chromosphere and corona 


(iMclntosh et all _ 

for the occurrence of ma gnetic reconnection processes (IZweibel fc YamaHa . 20091) 


20121) are some suggested mechanism s. There is ample evidenc e 


magnetic wave dynamics (INakariakov fc Verwichtd. 12005 ); and mass flows connect¬ 


ing chromospheric and coronal regions (|Peter fc Judge! . Il999l) . All these processes 
are under theoretical and observational study and advances on their description 
have been made so a discussion in terms of mutually exclusive mechanisms is be¬ 
side the point. It is important however to assess the importance of each mechanism. 

Wave energy transport and dissipation were regarded as relevant soon after the 
coronal heati ng problem aros e. Th e first proposed theo ries involved wave-based 
mechanisms. iBiermaniil ( 19461) and Schwarzschildl ( 19481) suggested that acoustic 
waves, generated in the convection zone, supplied the non-radiative energy needed 
to heat the chromosphere and corona in the form of m echanical energy trans¬ 
ported outwards (see review by Stein fc Leibacheil 19741) . Magnetic wave heating 
was proposed by Alfven ( 19471) soon after the ex i stence of magnetohydrodynamic 
(MHD) waves was postulated by himself ( Alfvenl . 1942 ). It was soon realised that 
the dissipation of magnetic waves could have an important role in coronal heat- 
i ng. because su rface Alfven waves ca n be transmitte d upwards in the atmosphere 
( Hollwed . 1978 ) and heat the corona ( Wentzell . 19741). These wa ves are indeed sub- 
i ect to mechanisms, such a s resonant absorption (|lonsoil 1978h and phase-mixing 


( Hevvaerts fc Priest . 19831) . that enhance the dissipation of energy by resistive or 


viscous processes, which would otherwise take too long in the high Reynolds num¬ 
ber co rona. Theoretical and numerical studies on the dissipation of hydro magnetic 
19791) and the resonance absorption of Alfven waves (lonson. 


waves 


(IWentzel. 


1978tlRae fc Roberts 


Poedts et al. 


lance aosorp ' 
1989Lll99d) 


were carried out in open and 


closed magnetic structures. Numerical experiments included the sideways or foot- 
point driving of the structures, t r ying to mimi c the b uffeting of the solar surface 


( Goedbloed fc HalberstadtI . 

observ ations were made by 


1994: Ofman et al. 


Hollweg fc Sterling 



The first relations t o loop 
to explain data by iGolubl 


( 19801) . However, proper instrumentation able to detect waves and reliably measure 
their properties was not available. 

This changed near the end of the last century, when routine observations of 
the solar atmosphere made with instruments onboard the Solar and Heliospheric 
Observatory (SoHO) and the Transition Region and Coronal Explorer (TRACE) 
led to a number of disc overies. Remarkable examples are the transverse oscil¬ 
lations in coronal loops ( Aschwanden et 1999t Nakariakov et al\ . Il999l) and 


their rapid time damping; the compressive waves in short and long coronal loops 


( Berghmans fc Glettel . 19991 : Nightingale et ai . 19991 : De Moortel et ai . 20001) : the 


quasi-periodic oscillations of the Doppler shift and i ntensity in emission lines in 


coronal loops (|Kliem et all 120021: IWang et ali . 120021) and their damping; or the 
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quasi-periodic disturbances in op en magnetic regions of the corona (jOfman et al 


1997 : DeForest &: Gurmanl. 1998). Extended reviews on thes e observations can b e 


found in Nakariakov fc Verwichte (l2005h . De Moorteil ( 2005h . Aschwanden 
and De Moortel fc NakariakovI (|2012li . 

The observed wave activity is interpreted in terms of standing and propagat¬ 
ing MHD waves. Nowadays, the existence of t heoretically predicted MHD wave 
types in the solar atmosphere is well established ( Robertsl l200(i IXschwandeil 2003t 
Nakariakov fc Verwichtel . 120051: iBallestei , 2000[l. This has lead to t he development 


of solar atmospheric seismology dUchida 1970l : Roberts et all 1984) , which enables 
the determination of physical parameters in the solar atmosphere from the compar¬ 
ison of observed and theoretical wave properties. The application of this method 
to the solar atmosphere has been successful, but the wave heating perspective was 
discontinued for more than a decade. 

A renewed interest has now grown on the wave heating question. This was 
motivated by the evidence for ubiquitous wave activity at different layers of the 
solar atmosphere, gathered by observations with instruments onboard e.g., the 
Hinode and Solar Dynamics Observatory (SDO) spacecraft or the Hi-C rocket 
experiment. These high resolution observations have demonstrated the presence 


(De Pontieu et all 2007): coronal magnetic loons (McIntosh et all 20 III) nromi- 

nencenlasmas (Okamoto et at. 2007: Arregui et a,l.. 20121) 

; X-ra-' 

i^iets (Cirtain et al.. 

2007 

): and extended regions of the corona (Tomczvk et al\ 

I 2 OO 7 I 

Tomczvk & McIntosh 


20091) . The observed disturbances display time and spatial damping possibly indi¬ 


cating in situ energy dissipation. In some cases, their energy content seems to be 
large enough to compensate estimated energy losses. These observations have re¬ 
placed coronal wave heating at the forefront of the discussion, leading to the need 
to reasse ss the significance of wave energy transport and dissipatio n to the heating 
problem ( Cargill fc de Moorteil l20lit Parnell fc De Moortel . 2012 ). 

We discuss the role of magnetic wave energy transport and dissipation processes 
in the heating of the solar corona, based on our current theoretical understanding 
and observational evidence. We focus on the following aspects: in Section [SJ we 
present a selection of wave activity observations by pointing out their relevance to 
the heating problem; Section[3]offers a description of theoretical models and physical 
processes that have been suggested to explain observations; particular emphasis is 
made in Section 0] on the series of physical processes that link wave energy content, 
propagation, damping and dissipation. Section[^contains a commentary on how the 
combination of data and theory, by means of forward modelling of observational 
signatures and seismic inversion of observed wave properties, should help advance 
the field in the future. Our conclusions are presented in Section [51 


2. Observed Wave Activity 


in solar coronal structures (see e.g., 

Billingsl 1959: Rosenberg. 1970l: Vernazza et a.l. 

1975: Tsubaki. 19771 Trottet et all 

I 979 I: Antonucci et al.. 1984 Aschwanden. 1987 

Deubner & Flecl^ 19891) . The detection of these oscillations was based on the mea- 


surement of the temporal and spatial variation of spectroscopic properties (intensity, 
width, and Doppler velocity) of coronal emission lines. Some of these observations 
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took advantage of the instrumentation onboard Skylab or Yohkoh spacecraft, but 
were mostly re stricted to time series analyses with little spatial information (see 
Tsnbakil . 1 19881 for a summary of early observations). 

This situation changed with the advent of the imaging and spectroscopic instru¬ 
ments onboard SoHO and TRACE spacecraft. Due to the temperature discrimina¬ 
tion and spatial resolution capabilities of their EUV and soft X-ray telescopes, ob¬ 
servations demonstrated the existence of wav e-like dynamics in the form of e.g., den¬ 


sity f l uctuations in plume s and coronal holes ( Ofman et al . , 1997 ; DeForest fc Gnrmanl . 


19981: Ofman et al. 2000l). propagating compressional disturbances in coronal loops 


2001 

De Moortel et al. 

19991 

Nakariakov et al. 




2002^ ol ), or transverse coronal loop oscillations (jAschwanden et al\. 


199911 . 


The existence of waves and oscillations in magnetic and plasma structures of the 
solar atmosphere is now beyond question. The spatial, temporal, and spectral reso¬ 
lution of imaging and spectroscopic instruments in current ground- and space-based 
observatories (SST, DST, SoHO, TRACE CoMP, Hinode, STEREO, SDO, HI-C, 
IRIS) have enabled us to directly image and measure motions associated to wave 
dynamics with increasing precision. Data show evidence for waves in increasingly 
finer plasma structures belonging to regions of the solar atmosphere with differ¬ 
ent physical properties, such as coronal loops; prominence plasmas; chromospheric 
spicules and mottles; coronal holes, plumes, etc. Let us focus on some particular 
examples and draw attention to the aspects that are relevant to wave heating. 

Transverse loop oscillations offered the first imaging evidence for stand ing waves 


in the corona associated with the periodic displacement o f these structures ( Aachwanden et al 


19991 : Nakariakov et all . Il999l: Aschwanden et al . , 20021 : Schriiver eA al . , 200211 . Iii 


these events wave damping can be directly observ ed and measured. This constitutes 


impo rtant information for seismology diagnostics (lArregui et al.l . l2007HPascoe et al 


2013f) and the reason for co nsidering damping mechanisms such as resonant absorp¬ 
tion ( Goossens et all l2006ll . However, TRACE data showed that these were occa¬ 
sional events and loops are hot irrespective of the presence or absence of lateral 
displacements. On the other hand, the Coronal Multi-channel Polarimeter (CoMP) 
observations indicate that coronal disturbances are present in extended regions of 
the corona (1.05-1.35 Rsun)- The measured Doppler velocity fluctuations, of the or¬ 
der of 0.3 km s“^, do not produce significant intensity variations and are interpreted 
as Alfven waves propagating along th e coronal magne t ic field, with measured phase 


speed s of the order of 1 to 4 Mm s ^ ( Tomczvk et al . , 2007 ; iTomczvk fc Mclntoshl . 


2 OO 9 II . Importantly, CoMP data show signatures of in situ wave damping in the form 
of a discrepancy in the outward to inward wave power, but energy estimates (~ 0.01 
W m“^) seem to fall below the amount required to heat the ambient plasma. The 
Solar Dynamics Observatory (SDO ) has shown us that Al fven waves are common in 


the transition regions and corona (IMcIntosh et aLl. l201lh and t hat motions visible 
in both regions share a common origin (jPe Pontien et al. . 2011 b Energy estimates 
vary from one region to another. The disturbances seem to be energetic enough 
to power the quiet Sun corona and coronal hole regions, with an estimated wave 
power of ~ 100 - 200 W m“^, but not the active region corona (~100 W m“^ in 
front of the requ ired 2,000 W m~^). Observati ons of transverse waves in the active 
region corona by Morton fc McLanghlinI ( 2013 b obtained by combining the highest 
available resolution instruments onboard Hi-C and SDO/AIA, confirm that wave 
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activity is all-pervasive in the EUV corona. However, this activity is concluded to 
be of low energy because of the relatively low displacement and velocity amplitudes 
(< 50 km and < 3 km s“^, respectively) in the measured waves. 

Chromospheric spicules are excellent candidates for mass and wave energy trans¬ 
port to the corona, since they provide a physical connection between both me¬ 
dia. They are present everywhere and at all times and show a complex dynamics 
with a combination of la t eral swaying, upward mass flo ws, and torsional motions 
(lOkamoto fc De PontieuL 20 lit De Pontieu et al. . 2012 tl . Recent observations by 
iMortonI ( 2014ll suggest that wave damping is present in these structures. The recent 
analyses by the IRIS team will clarify m a ny aspects about thei r field and plasma 
structuring and dynamics ( Pereira et all [2014 : Skoesrud et all . 2014). Propagat¬ 
ing and standing transverse oscillations are also present in chromospheric mottles 
belonging to strong photospheric field concentra tions. By measuring t he amplitude 
and phase speed variation along these structures, Kuridze et al. ( 2013ll suggest that 
their sudden change at a given height is an indication of the presence of non-linear 
wave processes. Photospheric field concentration s, such as bright p oints, were pro¬ 
posed as the source of torsional Al fven waves by Jess et al. ( 20091). The increasing _ 

evide nce for surface vortex motions ( Bonet et at. . l200^ Wedemever-Bdhm fc Rouppe van der Voort . 
2 OO 9 II strongly suggests this type of dynamics c an be important in the ge nera- 


gen 

tion and wave energy propagation to upper layers (iMathioudakis et al\ . l2013h . The 


chromosphere seems to be a vast reservoir of MHD wave energy with observations 
showing a complex combination of intensity and width v ariations of ch r omos pheric 
structures together with their transverse displacement ( Morton et a/.l . 2012 1. The 
evidence poi nts to the photosphe ric excitation of the dynamics observed at the 


upper layers (|Morton et a,l\ . 1201311 


Waves in polar coronal hole regions are gaining ample attention. The mea¬ 
surement of wave characteristics, such as the amplitude, periodicity, and phase 
speed provide plasma and held diagnostics and c onstraints on theoreti cal mod¬ 
els of coronal heating and solar wind acceleration ( Baneriee et all . l201lh . Recent 
analyses have enabled the measurement of the spatial variation of the non-thermal 
velocity for Alfven wa ves and the electron density along the propagation direction. 
iBaneriee et al. ( 2009l l hnd that the non-thermal velocity is inversely proportional 
to the quadratic root of the electron density, in agreem ent with what is pr edicted 
for undamped radially propagating linear Alfven waves. Hahn et all ( 2012 1 report 
that the widths of spectral lines, which are assumed to be proportional to the wave 
amplitude, decrease at relatively low heights. This is interpreted as an indication of 
wave damping. The energy dissipated in between 1 and 1.3 solar radii is estimated 
to account for up to 70% of the amount required to heat the polar coronal hole 
and accelerate the solar wind. Similar analyses have qu antified the energy carried 
out and dissipated b y Alfven waves in these structures ( Bemoorad fc Abbol . 12012 : 


Hahn fc SavinL 2013 ). By combin ing spectroscopic measurements with a magnetic 


field model, Hahn fc Savinl ( 2014^ were able to trace the variation of the non ther¬ 
mal velocity along the magnetic field above the limb of a quiet sun region. They 
find that the waves are dissipated over a region centred on the top of the loops. The 
position along the loop where the damping begins is strongly correlated with the 
length of the loop, implying that the damping mechanism depends on the global 
loop properties rather than on local co llisional dissipati o n. Tr ansverse waves seem 
to be present also in solar plumes and iThurgood et 'al. ( 2014 1 have measured the 
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occurrence rate of wave properties, such as periods and displacement and velocity 
amplitudes in great detail. 

The relationship between wave signatures at different heights of the solar atmo¬ 
sphere is reported in many s tudies (see e.g., O’Shea et al . 2002 : Brvnildsen et oil 
2003; Marsh fc Walshl . 200till . They provide information about the phase speed of 
propagating waves fro m the time delay of signatures de tected at different heights in 
the solar atmosphere ( Abramov-Maximov et al. . 201lll . The high spatial and tem¬ 
poral resolution of SDO/AIA is enabling additional analyses of the spatial distri¬ 
bution of the frequency, time seri es analyses, and c orrelations betwe e n the signals 
at different atmospheric heights ( Reznikova et 120121: iFreii of. all l20l4l . How¬ 
ever, obtaining a causal connection between wave dynamics at different regions 
of the atmosphere is still problematic. Understanding how the waves are gener¬ 
ated and behave as a funct ion of the line forma t ion te mperature and the magnetic 
field structure is essential ( Mariska fc Mugl achl 2010 ) and observations should be 
complemented wit h numerical modelling (iKhomenko fc Calvo Santamaria . 2013 : 


Felipe et aH 120141 . 


3. Theoretical Models and Physical Processes 

The observed wave activity is interpreted as the manifestation of the presence of 
MHD waves. There are three basic MHD wave-types in a uniform plasma of infinite 
extent and they have properties that make them suitable for wave energy transport 
and eventual plasma heating. Fast waves transport energy across magnetic surfaces. 
They resonantly couple to Alfven and slow waves in non-uniform plasmas, which 
leads to the transfer of energy from large to small spatial scales. Slow waves are 
basic in the generation of wave dynamics in the lower solar atmosphere, where 
the plasma still plays the relevant role in front of the magnetic forces. They are 
crucial in non-linear wave transformation processes leading to the development of 
shocks and their eventual dissipation. Alfven waves can connect and carry energy 
between remote regions of the atmosphere. They propagate to large distances along 
the magnetic field lines. Their energy is not easily dissipated, but their dynamics 
is characterised by large cross-field gradients in non-uniform plasmas, which leads 
to the enhancement of viscous and resistive dissipative processes. 

None of t he observed waves is expected to happen with the propert ies described 
in textbooks ( Goossens . 2003t Goedbloed fc Poedtsl . 2004 Priest! 2014 ). Pure MHD 
waves do not exist in the solar atmosphere. The classic description in terms of fast, 
slow, and Alfven waves treated separately is not accurate enough for the descrip¬ 
tion of the observed dynamics. The highly non-uniform and dynamic nature of the 
plasma in the solar atmosphere leads to complex wave properties. We identify three 
levels of complexity in the theoretical and numerical modelling of MHD waves in the 
solar atmosphere: (1) the analysis of linear waves in simple magnetic configurations 
enable us to study basic properties of wave trapping and propagation and to un¬ 
derstand in detail the workings of individual physical mechanisms. They offer local 
information on plasma and field properties from seismology analyses, but observa¬ 
tional signatures with little potential to be compared to real data; (2) the analysis 
of non-linear waves in structured and dynamic plasmas enables us to study wave 
transformation processes, wave-flow interactions, and the modelling of relevant pro¬ 
cesses such as shocks waves, non-linear Alfven waves, and Alfven wave turbulence. 
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They produce useful synthetic data to extract conclusions from the comparison 
to observations; (3) global numerical models use prescribed or observed boundaries 
and drivers as an input on large-scale numerical simulations. They make possible to 
include almost any physical ingredient, to perform global seismology, and to obtain 
useful observational signatures and predictions. 

The three approaches offer valuable information and are highly complementary. 
Results from simple models cannot be directly used to extract solid conclusions, but 
offer detailed understanding about physical processes that might occur under real 
Sun conditions. Their benefit is that the mechanism(s) of interest can be studied in 
isolation and under controlled conditions. On the other extreme, sophisticated sim¬ 
ulations enable the direct comparison to observed data, but detailed understanding 
about the physical processes is lost when many of them are allowed to happen si¬ 
multaneously in complicated simulations. The level of detailed understanding will 
decrease as we increase applicability and vice versa. 

Resonant damping of oscillations is a well studied mechanism for energy trans¬ 
port from large scale transverse motions to localised small scale m otions and is due 
to th e inhomogeneity of the medium in the cross-field direction Isee lGoossens et aT 


2006 , for a review). The mechanism has been mai nly studied in flux tube models 


( R.uderman fc Erdelvil l200^ Goossens et a,l\ . I 2 OIIII . but a cylin drically symmetri c 


flux tube tube is not necessary for the process to operate (see e.g., Rial et al . 2013h . 


Resonant absorption is robust in front of model complicatio ns, such as non-linear 

evolution (Terradas et al. , 2008q ) or multi strand ed structure (iTerradas et all 2008 _ 

Pascoe et all . 2011 1. and works for both stand ing ( Goossens et al . 2002t Ruderman fc Roberts . 


2002) and propagating ( Pascoe et [ 2 ^ waves, producing time and spatial 
damping of the wave amplitude, respectively. This is an ideal physical process , 
indep endent of resistivity in the limit of high Reynolds numbers ( Poedts fc Kerneil 
I 1991 I) . Once in the resonant layer, energy can further scale down to smaller spa¬ 


tial s cales by the process of phase mixing ( Hevvaerts fc Priest . 19831 : iHood etd. 


20021) . The stronger the cross-field i nhomogeneity, t h e quic ker is the damping, but 


the slower the small scale creation (ITerradas et a, 11 1200611 . When a given length- 


scale is reached, viscous and/or resistive dissipation processes turn on, thus enabling 
the heating of the plasma. 

To date, two arguments support the mechanism of resonant absorption being 
operative in the observed waves. First, theory predicts damping time and spatial 
scales that are compatible with those observed. Second, resonant absorption is a 
frequency selective process ( Terradas et al ],[2H3), with low-frequency waves being 
less damped in front of high frequency w aves. A compari s on be tween the outward 
to inward power ratio measur ements by Tomczvk et aTI ( 20071) using CoMP and 
the theoretical predictions by IVerth etajJ ~i 2010 l) strongly supports the idea that 
resonant absorption might be producing the observed in situ ener gy loss. _ 


Another mechanism studied in detail is Alfven wave turbulence (Ivan Ballegooiien et al 


2011 . 20i 4 Asgari-Targhi et ^ . 120131) . According to the models, photospheric foot- 


point motions transported along loops get amplified and reflected at the transition 
region boundary, thus producing a pattern of counter-propagating perturbations. 
Their complex interactions lead to the creation of small scales and dissipation. 
Recent numerical experiments show that, if that is the case, the obtained dissipa¬ 
tion rates are able to reproduce chromospheric and coronal heating requirements. 
The heating rate scales with the magnetic field strength and with the loop length. 
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The modelling of the lower atmosphere is also an important factor. Realistic lower 
atmosphere modelling favours AC heating in front of DC dissipation. 


Ob servations are now starting to offer support in favour of this process. iDe Moortel et al. 
( 2014h analyse CoMP Doppler shifts oscillations and their frequency distribution 
at both sides of a large trans-equatorial loop system and find that an excess of high 
frequency power is present at the apex. This is interpreted as being due to low and 
mid frequency wave energy cascading down because of Alfven wave turbulence. A 
recent comparison between observed non-thermal velocities with predictions from 
theoretical models in coronal loops bv lAsgari-Targhi et al. (2014) leads to the con¬ 
clusion that Alfven wave turbulence is indeed a strong candidate for explaining how 
the observed loops are heated. 

Significant efforts are being done in the design and ap plication of numeri¬ 
cal models, such as the Alfven Wave Solar Model (AWSoM, ISokolov et al. 2013 : 
van der Holst et al. 2014ll . that address the coronal heating problem and the solar 
wind acceleration from a global perspective. These models incorporate Alfven wave 
turbulence. The result is a remarkable capacity to reproduce the observed EUV 
emission and produce solar wind predictions. 

Physical processes discussed above are mostly restricted to MHD models. How¬ 
ever, the MHD approximation is possibly part of the difficulty to advance the 
wave heating question, since the actual heating processes must occur at kinetic 
scales and hybrid or fully kinetic analyses need to be further developed. The ideal 
MHD model should fail when appli ed to the lower solar a t mosp here where par¬ 
tial ionisation effects are important ( Khomenko fc ColladosL 2012 ) and effects like 
ambipola r diffusion, collisions and non-magnetisation become non-negligible. For 
instance. IVranies fc Poedt^ ( 20101) have pointed out that, because of the presence 
of non-magnetised ions in the photosphere, Alfven waves cannot be efhciently gen¬ 
erated in, nor travel through, this region. More importantly, the wave energy flux 
through the photosphere becomes orders of magnitude smaller, compared to the 
ideal case, wh en the effects of part ial ionisation and colli sions are consistently taken 
into account (Vranies et al . 2008h . On the other hand, ISong fc Vasvliunirs ( 201lh 
and ITu fc Song! (120131) have studied Alfven wave propagation and heating through 
plasma-neutral collisions and find that they can generate sufficient heat, with most 
of the heat deposited as required at lower altitudes. To clarify these aspects is cru- 
cial in view of the recen t observations of the chromospheric dynamics with IRIS 
2014 . 


(De Pontieii et al 


4. Wave Energy Content, Damping and Dissipation 

Our estimates for the energy content on the observed wav es are based on simple 
expressions for bulk Alfven waves in homogeneous plasmas. iGoossens et ali (120131) 
have shown that the energy flux computed with the well known expression for bulk 
Alfven waves could ove r estim ate the real flux by a factor in the range 10 to 50. 
Van Doorsselaere et al. ( 20l4 provide approximations to the energy propagated 
by kink modes in an ensemble of flux tubes finding a correction factor for the 
energy in kink waves, compared to the bulk Alfven waves, in terms of the density 
filling factor. 

The full energy balance will be given by Poynting theorem, which relates the 
amount of electromagnetic energy variation; the fraction of that energy that prop- 
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agates with the wave; and the remaining energy that is dissipated. The assessment 
of these quantities requires the measurement of the spatial and temporal variation 
of both velocity and magnetic field perturbations. This being a challenge, a first 
step is the detailed analysis of the flow of energy and its spatial distribution in the 
waveguide. For reso nance absorption in a flux tube model this has been done by 
Arregui et al. (201 ill . in the prominence context. The results show that ideal damp¬ 


ing is produced by the jump of the radial component of the Poynting vector, which 
produces the energy flow into the dissipative layer from both sides of the resonance. 
Once in the resonant layer, energy flows along the field and resistivity, wave fields, 
and the created currents determine the energy dissipation and its spatial distribu¬ 
tion. Quantification and comparison to observed field-aligned brightness variations 
remain to be done. 

In numerical experiments, all magnetic and plasma perturbation are available. 
The kinetic, magnetic, and total energy and their spatial and temporal evolution 
can readily be analysed. This is not the case with observations, where the energy 
contained in t he magnetic field is hidden a nd has to be inferred. The numerical 
simulations bv IPe Moortel fc Pascod ( 2012 b for transverse wave propagation in a 
multi-stranded waveguide model, have shown that line of sight integration effects 
can be important when evaluating the energy content of a wave from only the 
measured kinetic energy. Also, that a relevant amount of energy that is present in 
the form of magnetic perturbations in the simulation would be out of observational 
reach. The “visible” energy would only be in the 3% to 10% range of the total 
energy. 

As first noted by iLee fc RobertsI (|l980li , wave damping is not a synonym for 
wave dissipation. The two processes might operate at different time and spatial 
scales and e.g., viscosity should not immediately dissipate the enhanced local os¬ 
cillations. The matching of the two time-scales is important. Consider for instance 
resonant damping, phase mixing, and resistive dissipation. A calculation by J. Ter- 
radas (personal communication) shows that, under coronal conditions (magnetic 
Reynolds number of the order of 10^^) the observational consequence, the attenua¬ 
tion of the motion, might give little information on the concealed physical process, 
the dissipation of the wave energy. For instance, for a typical coronal loop subject to 
resonant damping because of a transverse inhomogeneity length scale of l/R = 0.1 
(in units of the tube radius i?), the damping time over per period is t^/P ^ 13. The 
time at which dissipation becomes important, once phase mixing has sufficiently 
decreased the characteristic transverse spatial scales, is r^afP ^ 170. The stronger 
the non-uniformity, the faster is the damping, t^/P ~ 3 for l/R = 0.5, but the 
slower the small scale creation due to phase mixing. We need therefore the wait 
more time for resistivity to become important, Tra/P ^ 500. This means that en¬ 
ergy dissipation would only be efficient at very large times, and there would be no 
heating during the observed oscillations. 

The numerical experiment on resonant absorption bv lOfman et al. (1998) anal¬ 
ysed the plasma response and the creation of small-scale structure because of the 
non-linear response of the density, but the assumed energ y balance was n o t entire h 


self-consistent. More recent numerical experiments by iTerradas et al\ (|2008 


arely 


have analysed the first two stages (damping and phase mixing), in complex mod¬ 
els including the non-linearity, but no heating was computed. A detailed numerical 
study of the full process of wave damping, phase mixing, resistive or viscous dissipa- 
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tion, and plasma heating, including the forward modelling of the synthetic plasma 
response has not yet been undertaken. Such an experiment would prove/disprove 
the viability of heating by resonant absorption and phase mixing, by producing 
spatial and temporal heating profiles to be compared to observations. 

5. Confronting Theory and Observations 

Most of the efforts in the area have focused on the design and use of better in¬ 
strumentation for data acquisition and analysis and the improvement of theoretical 
and numerical modelling. The two developments being important, further advance 
will only come from the development of informative tools for the combination of 
data and theory to assess the physical conditions and processes that operate in 
the solar atmosphere. The combination of observed data and theoretical results to 
infer physical properties and mechanisms is not easy task. One needs to solve two 
problems simultaneously. In the forward problem, we prescribe theoretical mod¬ 
els and parameters (the causes) and analyse the theoretical wave properties (the 
consequences). In the inverse problem, we try to infer the causes (the unknown 
physical conditions/mechanisms) from the consequences (the observed wave prop¬ 
erties) . The discovery of wave dynamics pervading the solar atmosphere has led to 
considerable efforts in both forward and inverse modelling. 


(a) Forward Modelling of Observational Signatures 

The creation of synthetic imaging and spectral observational signatures for di¬ 
rect comparison to data is starting to be used as a means to extract conclusions 


about the goodness of our wave mode l s Isee e.g..lSh elvag et a/.l. 120061: iTarova n et al^ 
20071 : IPe Moortel fc Bradshaw . 2008; Tarovan fc ErdelviT 2009^ ' Owen et oK 20091) . 


Forw ard modelling is essential to solve wave type identification issues. iVan Doorsselaere et a,I 
( 2008ah raised the concern as to whether the observed transverse motions should 
be interpreted as fas t kink waves or Alfven waves. The issue was clarified by 
Goossens et a,l. ( 2009ll . who pointed that the observed dynamics is the result of 
Alfven waves modified by the radial structuring of the plasma density. Their nature 
being highly Alfvenic, these authors refer to them as surface Alfven waves. Kink 
and torsion al Alfven waves in ma gnetic flux tubes produce distinct observational 


signatures (jGoossens et a/.l . l20l4) . Besides the naming, which is of secondary im¬ 


portance, it is crucial to consider an adequate description of the structures in which 
the waves propagate, in particular whe n it comes to se i smolo gy and energy budget 
calculations using the observed waves ( Goossens et oL . 20I2ll . Observed tilted fea- 
tures in spectrograms ha ve been interpreted as sign atures of torsional Alfv en waves 
( De Pontieu et al . 2012 ). The forward modelling bv iGoossens et al. ( 20141) suggests 
that tilted features may also arise under the kink wave interpretation, depending 
among other factors on the line of sight integration. 

Another applications of forward modelling techniques include the fitting of wave 
properties, such as wave periods and amplitudes, frorn the co mparison of alternative 
data realisations and observed data ( McIntosh et al. . 201 ih : the assessment of the 
wave or flow nature of observed q uasi-periodic disturbances from Doppler shifted 
spect ra ( De Pontieu fc Mclntoshl . l201(t iTian et d. . 2011 . 2012 : Verwichte et a,l\ . 
2010l:lKiddie et ci.l.l2012l) : or the analysis of the nature of shock distributions along 
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coronal loops as bein g due to either non-linear Alfven wave or nano flare heating 
( Antolin et al\ . 2008ll . 

The potential of forward modelling should still be further exploited to compare 
alternative wave models producing distinct observational consequences and alter¬ 
native heating mechanisms. Different heating models are known to produce distinct 
heating profiles along waveguides. For instance, viscous heating is determined by 
the velocity perturbation profile along waveguides, leading to apex heating, while 
resistive dissipation is determined by magneti c field perturbations, thus produc¬ 
ing foot-point heating ( Van Doorsselaer e et a.l\. 2007tl. The recent developme nt of 
forward mode l ling t echniques to waves by Antolin fc Van Doorsselaei^ ( 201,'lh and 


Antolin et al\ (1201411 show great promise. 


(6) Seismology Inversion from Observed Wave Properties 


In parallel to theoretical and observational advances, the field of solar atmo¬ 
spheric seismology has emerged. The term refers to the study of the physical condi¬ 
tions in solar atmospheric magnetic and plasma s tructures from the analysis of the 
observ ed wave properties. It was first s uggested by Uchida ( I97C 1 and iRoberts et al. 


(1984), in the coronal context, and bv iRoberts fc Joardei ( 1994 1 in the prominence 


context. The aim is to increase our knowledge about the complicated structure and 
dynamics of the solar atmosphere and is based on the fact that the properties of the 
observed oscillations are determined by the plasma and magnetic field pr operties. 


2000 

iDe Moortel. 

2 OO 5 I 

Nakariakov & Verwichtel. 20051 Baneriee et a/.L 120071: lArregui. 

2012 

iDe Moortel & Nakariakov, 

2012 

, for reviews). This was made possible by the 


increase in the quantity and quality of wave activity observations and the refinement 
of theoretical MHD wave models. 

Local seismology, based on the use of simplified models for the magnetic and 
plasma structuring of tube-like waveguides and observations of the period and 
damping of standing and prop agating waves, have enabled us to obtain information 


on the magnetic field strength (INakariakov fc O fmanl.l2001l:lyan Doorsse laere et al. 
l2008fJl , the coronal density scale heigh t (jAndries et ah . l2005l 12009 ), the mag¬ 


netic field expa nsion (IVerth et all 1200811 . the Alfven speed (lArregui et all 12007 


Goossenj. 20081 or the longitudinal and cross-field magnetic field and density struc¬ 


turin g ( Verwichte et al . 2006t Verth et al. . 2011 : Arregui et al\ . 20136t Arregui fc Asensio Ramos 

I 2 OI 4 II . This information is crucial to evaluate the time and spatial scales for wave 
damping and dissipation processes (Section 0]). The simplicity of the models con- 
sid ered so far imposes a limited applicability to real Sun conditions. As pointed out 
bv IPe Moortel fc Pascod ( 2009f l. current inversion results need to be re-examined 
in front of the outcome from more involved numerical models, in order to determine 
their accuracy. A s the number of obser ved events incr eases, statistical analyses - 


both frequentists (IVerwichte et al.l . l2013ri and Bayesian (lAsensio Ramos fc Arregui 


2 OI 3 II - offer valuable wide-ranging information. 


Seismology diagnostics will play a key role when it comes to quantify the energy 
carried out and dissipated by MHD waves, but the area needs to mature by consid¬ 
ering models more akin to the structure and dynamics of solar plasmas. Aiming at 
the inference of magnetic field and plasma properties in spatially extended regions 
of the atmosphere is another important topic. The available maps with measured 
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wave characteristics by e.g., Tomczvk fc McIntosh ( 20091) offer a way to perform 
glob al seismology , but few global inv ersion techniques have been devised so far (see 
e.sf.. lBallaiL 2007 : Kwon et ai . 2013). 


(c) Bayesian Analysis 


Because any inversion process has to be done under conditions in which infor¬ 
mation is incomplete and uncertain, probabilistic inference is the next step. The 
solution to a probabilistic inversion proble m is given in the form of statements 
about the plausibility of parameters/models ( Javnesl . 2003 ). Bayesian analysis then 
enable us to quantify the degree of belief on these statements, by measuring to 
which amount their are supp orted by information on observed data. This is done 
by application of Bayes’ rule ( Baves fc Pricel 1763 ). a mathematical theorem that 
teach us how to combine prior knowledge with the likelihood of obtaining a data 
realisation as a function of the unknown parameters/models to compute the so- 
called posterior, which accounts for what can be said about a parameter or model, 
conditional on data. 

The concept can be applied to the problems of parameter inference and model 
comparison. In parameter inference the posterior is computed for different combina¬ 
tions of parameters and then one marginalises the full posterior to obtain informa¬ 
tion about the parameter of interest. In model comparison, the ratio of posteriors 
for alternative models is computed to assess which one better explains observed 
data. Bayesian analysis is producing successful results in the physical sciences (see 
von Toussaintl . l20Ill) and other areas of spac e science and astrophys i cs re- 


e.g _ _ _ . _ _ _ 

search such as cosmology or exoplanet detection ( Loredol . Il992l : I Gregor 3 , 1200,4 
Trottal . 120081 ) . Applications to sol ar physics ar e scarc e, but the few applications 


to the analysis of sol ar oscillations (Marsh et al . 2008h . the comparison of heating 


profiles along loops (lAdamakis et ali . 120101) or to inference and model compar¬ 


ison problems in sola r atmospheric seismology ( Arregui &: Asensio Ramo^ . l201ll 
Arregui et all l2013al ) show great promise. 

The main advantages of the Bayesian methodology are that inference is made 
using self-consistently all the available information in prior knowledge, observed 
data, and model, additionally obtaining a correct propagation of uncertainty. In 
model comparison, the method enables to assess quantitatively which model among 
competing alternative explanations better explains observed data. The development 
and application of Bayesian techniques can shed light to the problems of wave type 
identification, wave vs. flow interpretation, and the comparison between alternative 
wave damping and heating scenarios. 


6 . Conclusions 

Physical explanations based on waves were at the forefront of the coronal heating 
discussion since the problem came to existence. At that time, evidence about the 
presence of waves and oscillations was absent and the role magnetic fields could 
play was largely unknown. The solar atmosphere is now one of the best observed 
and studied astrophysical systems and the presence of waves and oscillations is 
beyond question. Waves and oscillations are found and analysed in structures with 
very different physical conditions. Multi-wavelength and multi-scale imaging and 
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spectroscopic observations enable us to measure their properties with increasing 
precision. 

Wave heating theories remain plausible, but most current models are still too 
simple to be applicable to the real Sun. Future models should consider the highly 
structured and dynamic nature of the coronal plasma. The observed waves act as 
energy carriers and can transport energy to small scales, where the heating pro¬ 
cesses occur concealed from observational scrutiny. The theoretical analysis of such 
kinetic processes should be pursued to derive their large scale observational conse¬ 
quences. Similarly, the consequences of partial ionisation and collisions on the wave 
properties and their energy dissipation need to be quantified. These observational 
consequences need to be fully developed, by forward modelling of synthetic data and 
comparison to observations. As alternative explanations for the same phenomena 
arise, there is a need to devise model comparison tools to assess the performance of 
alternative wave energy transport and dissipation mechanisms in explaining data. 
Similar methods need to be formulated to confront wave based mechanisms with 
another explanations, such as mass flows or nanoflare heating. 

Wave energy seems to fill the solar atmosphere. Enough energy is available at 
the photospheric level and part of it is transmitted above. We do not know exactly 
how and in which amount. Future observations should concentrate on tracking the 
flow of energy across different regions of the atmosphere. The amount of energy in 
magnetic perturbations is hidden. Obtaining estimates of the energy in small scales 
and dealing with line of sight integration effects represent a challenge for data 
analysis. The solution is to combine the appropriate modelling with the application 
of forward modelling and inference tools. 

Advances in both theory and observations have placed again waves at the fore¬ 
front of the discussion. The design and application of tools for comparison between 
theory and observations is essential. Forward modelling is in its early babyhood and 
seismology is in its infancy. The development of a self-consistent methodology to 
combine information from theory and data, such as Bayesian analysis, is still await¬ 
ing to be developed and applied to this area. This would enable us to establish what 
can be plausibly said about the physical conditions and processes operating in the 
observed wave dynamics, by using all the available information. 

This work was supported by a Ramon y Cajal Fellowship and project AYA2011-22846 
from the Spanish Ministry of Economy and Competitiveness (MINECO). I am grateful to 
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